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Abstract 

The transition from phase chaos to defect chaos in the complex Ginzburg-Landau equation 
(CGLE) is related to saddle-node bifurcations of modulated amplitude waves (MAWs). First, 
the spatial period P of MAWs is shown to be limited by a maximum P$n which depends on 
the CGLE coefficients; MAW-like structures with period larger than P$n evolve to defects. 
Second, slowly evolving near-MAWs with average phase gradients v ~ and various periods 
occur naturally in phase chaotic states of the CGLE. As a measure for these periods, we study 
the distributions of spacings p between neighboring peaks of the phase gradient. A systematic 
comparison of p and P$n a s a function of coefficients of the CGLE shows that defects are 
generated at locations where p becomes larger than Psn- In other words, MAWs with period 
Psn represent "critical nuclei" for the formation of defects in phase chaos and may trigger the 
transition to defect chaos. Since rare events where p becomes sufficiently large to lead to defect 
formation may only occur after a long transient, the coefficients where the transition to defect 
chaos seems to occur depend on system size and integration time. We conjecture that in the 
regime where the maximum period Psn has diverged, phase chaos persists in the thermodynamic 
limit. 

PACS: 05.45. Jn 03.40.Kf; 05.45.-a; 
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1. Introduction 

The transition from phase to defect chaos for the one dimensional complex Ginzburg-Landau 
equation (CGLE) was recently related to the bifurcation properties of a family of coherent 
structures called modulated amplitude waves (MAWs) M. In this paper the relationship between 
MAWs and large scale chaos is studied in detail, providing a comprehensive description of various 
aspects of the CGLE chaotic dynamics. 

When a spatially extended system is driven sufficiently far away from equilibrium, patterns 
can eventually form |2|||. In many cases these patterns show an erratic behavior in space 
and time: such behavior is commonly referred to as spatiotemporal chaos 0-f|- Examples of 
extended systems displaying such chaotic dynamics in one spatial dimension include: heated 
wire convection H, printers instability and film drag experiments [pj, eutectic growth [||, binary 
convection ||, sidewall convection fllQfl , the far field of spiral waves in the Belousov-Zhabotinsky 
reaction [11|, the Taylor-Dean system |l2]| , hydrothermal [ 13 1 and internal 14] waves excited by 



the Marangoni effect and the oscillatory instability of a Rayleigh-Benard convection pattern []15| . 

Near the pattern forming threshold, the dynamics of such systems can often be described 
by so-called amplitude equations. When the pattern forming bifurcation from the homogeneous 
state is a forward Hopf bifurcation, the appropriate amplitude equation is the CGLE ||U, which 
in one spatial dimension reads as: 

d t A = A + (1 + ic x )dlA - (1 - ic 3 )\A\ 2 A , (1) 

where c\ and C3 are real coefficients and the field A = A(x, t) has complex values. 

For different choices of the coefficients numerical investigations of the CGLE have revealed the 



existence of various steady and spatiotemporally chaotic states [|]-|5,15-26|. Many of these states 
appear to consist of individual structures with well defined propagation and interaction proper- 
ties. It is thus tempting to use these structures as building blocks for a better understanding of 
spatiotemporal chaos. In this paper we will essentially follow such an approach. 

As a function of the coefficients c\ and C3, the CGLE (|]) can exhibit two qualitatively different 
spatiotemporal chaotic states known as phase chaos (when the modulus |A| is at any time 
bounded away from zero) and defect chaos (when \A\ can vanish leading to phase singularities). 
It is under dispute whether the transition from phase to defect chaos is sharp or not, and if a 
pure phase-chaotic, (i.e., defect-free) state can persist in the thermodynamic limit [21,27]. We 



will address these issues by suggesting a mechanism for the formation of defects related to the 
range of existence of MAWs. 

The main points of this paper are outlined in the following and illustrated in Figs. ||J2|. (i) 
Our investigation starts with the study of MAWs, which are uniformly propagating, spatially 
periodic solutions of the CGLE. These MAWs are parameterized by the average phase gradient 
v and their spatial period P. Our study is confined to the case v = for reasons specified below. 
Spatial profiles and the stable propagation of a particular MAW are presented in Fig. |l|a-c. 
Isolated MAW structures consisting of just one spatial period P play an important role in defect 
formation. In particular, for fixed CGLE coefficients the range of existence of coherent MAWs is 
limited by a saddle-node (SN) bifurcation which occurs when P reaches a maximal period Psn- 
(ii)If the MAWs are driven into conditions with P > Psn a dynamical instability occurs leading 
to the formation of defects (Fig. |l]d). (Hi) Slowly evolving structures reminiscent of MAWs 
( "near-MAWs" ) are observed in the phase chaotic regime (Fig. |l|e,f). In order to characterize 
such states, we have examined the distribution D(p) of spacings p between neighboring peaks of 
the phase-gradient profile. In particular for sufficiently long spacing p, the observed phase chaos 
structures are often very similar to a single period of a coherent MAW (Fig. |]f). (iv) When a 
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Figure 1. Summary of our main results which constitute a picture for the formation of defects 
from phase chaotic states. (a,b) Example of a coherent structure: phase gradient and modulus 
of a period P = 30 MAW at c\ = 0.6, C3 = 2. (c) Space time plot showing the stable propagation of 
the MAW from (a,b) in a small system of size P with periodic boundary conditions. Subsequent 
space time plots also show the phase gradient encoded in gray-scale (minima appear dark, 
maxima bright), (d) The same MAW as initial condition creates defects at c\ = 0.7, C3 = 2 
where P > Psn = 26.8. Black bars above the x-axis denote the size of Psn specific to the 
parameters of the panel. (e,g,h) Large scale chaos at c\ = 0.63, C3 = 2, L = 512. (e) Snapshot 
of the phase gradient profile with individual inter-peak spacings p. (g) Space time evolution 
of phase chaos and (h) distribution D(p) showing p <C Psn and no defects. A transient of 
t~10 is not shown. (f,i,j) Large scale chaos at c\ = 0.65, C3 = 2,L = 512. (f) Snapshot of the 
phase gradient profile t = 120 before the first defect forms and the MAW (dotted, P = Psn) 
overlayed onto the long structure, (i) Transient phase chaos with a fast and long structure 
traveling through the system which eventually nucleates defect chaos at t = 400, x = 360 (a 
transient of t ~ 10 is not shown). A snapshot of this structure was shown in (f). (j) The tail 
of the distribution of p reaches p > Psn due to the long structure; this leads to the break down 
of phase chaos. The distribution D(p) shown in (h) is also reported (dashed line). From the 
comparison of the two it is evident that the distributions do not modify dramatically when c\ 
is increased, while Psn decreases noticeably. 
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Figure 2. Phase diagram of the CGLE, showing the Benjamin- Feir-Newell curve (thin dot- 
dashed) where the transition from stable homogeneous oscillations to phase chaos takes place. 
The curves L\ (long dashed), L2 (thin dashed) and L3 (dashed) as obtained in [17, 19f| separate 
the various chaotic states. The filled circles correspond to our estimates of the L\ and L3 
transitions based on direct simulations of the CGLE along the 17 cuts in coefficient space that 
we studied. The open circles correspond to the location in coefficient space where the maximal 
inter-peak spacing p max is equal to the maximal MAW period Psn- Only small discrepancies 
between these two can be observed. Finally the full curve shows the Psn - * 00 limit which we 
conjecture to be a lower boundary for the transition from phase to defect chaos. 



phase chaotic state displays spacings p larger than Psn, phase chaos breaks down and defects 
are formed (e.g. at t = 400, x = 360 in Fig. P). Thus, the MAW with P = Psn may be viewed 
as a "critical nucleus" for the creation of defects. In phase chaos defect formation is similar to 
the dynamical process by which isolated MAW structures generate defects (Fig. |T|d). Therefore 
purely phase chaotic states are those for which p remains bounded below Psn (Fig. ||g), while 
defect chaos can occur when p becomes larger than Psn (Fig. [l|i). (v) A more detailed study of 
the probability distribution of the p's shows that for large p the probability decays exponentially 
(Fig. |l]h,j). As long as Psn has a finite value, we expect that, possibly after a very long transient 
time, defects will be generated, (vi) However, in a finite domain of the phase chaotic region, 
MAWs of arbitrarily large P exist: we expect that in this region, even in the thermodynamic 
limit, phase chaos will persist. Fig. [2] shows the main quadrant of the CGLE coefficient space. 
The region of persistent phase chaos is bounded by the Benjamin- Feir-Newell curve (thin dot- 
dashed) and the curve along which Psn ^ 00 (full curve in Fig. |2|). 

The outline of this paper is as follows: Section || is devoted to the study of the coherent 
MAW structures. In section 2.2 we study the bifurcation diagram of the MAWs, starting from 



the homogeneous oscillation. In section |2.3| the incoherent dynamics of near-MAW structures 
is presented. We show that for p>Psn, i-e., beyond the saddle-node bifurcation, near-MAWs 
evolve to defects. To illustrate the origin of the saddle-node bifurcations in section 2A we 
compare bifurcation diagrams of coherent structures for different phase gradient expansions of 
the CGLE. For the lowest order expansion (known as the Kuramoto-Sivashinsky equation Q) 
the saddle- node bifurcation is absent while it is captured by expansions of higher order. This 
explains why the divergence of the phase gradient was exclusively observed in simulations [20| 
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of higher order expansions. In section |3] we study various aspects of spatiotemporal chaos in the 
CGLE, and relate the observed continuous (L±) and discontinuous (L3) transitions (see Fig. |2|) 
to properties of the MAWs. The transition to defect chaos takes place when near-MAWs with 
periods larger than Psn occur in a phase chaotic state. In section [O] the typical values of p in 
the phase chaotic regime are related to the competition of two instabilities of the MAWs, and it 
is possible to give a good estimate for the numerically measured transition from phase to defect 
chaos from these considerations. A discussion of the presented results and some final remarks 
are reported in section ||. 



2. Modulated Amplitude Waves 

In this section we study the main properties of modulated amplitude waves (MAWs) 



First, in section 2.1 the coherent structure framework that we use to describe the MAWs is 



introduced. The bifurcation diagram of MAWs is explored in section 2.2, with a particular 



focus on the saddle-node bifurcations that limit the range of existence of MAWs. In section 



2.3 we study the nonlinear evolution of near-MAWs that are "pushed" beyond their saddle- 



node bifurcation and show that this leads to the formation of defects. Finally, in section \IA 
a bifurcation analysis of MAW-like coherent structures is performed in various phase equations 
that have been proposed as approximated models for the phase chaotic dynamics of the CGLE, 
and we show that only higher order phase equations reproduce the saddle- node bifurcation. 
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Figure 3. Examples of ODE solutions and corresponding amplitude and phase gradient profiles 
of MAWs. (a) Homoclinic orbit for ci=0.55 and 03 = 2; (b,c) corresponding profiles, (d) Limit 
cycle for c\ = 0.60, C3 = 2 and P = 30; (e,f) corresponding profiles. Dots in (a,d) denote the 
unstable fixed point (1,0,0) from which these orbits emerged. 
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2.1. Coherent Structures 

Coherent structures in the CGLE are uniformly propagating structures of the form 

A(x, t) = a(x - vty^-^e^ 1 , (2) 

where a and <p are real- valued functions of z:=x — vt. Coherent structures have been studied 
extensively [23-23] and play an important role in various regimes of the CGLE [ 16, 22 - 25| . 



The restriction to uniformly propagating structures reduces the CGLE to a set of three coupled 
ordinary differential equations (ODEs) |28|] . These ODEs are readily found by substitution of 
Ansatz (|2|) into the CGLE (|l|) and read as: 

a z = b 

b z = i^a. — 7 _1 [(1 — c\uj)a + v(b + citpa) — (1 — cics)a 3 } (3) 
-02 = —2hp/a + 7 _1 [ci + uj + v(cib/a - -0) - (ci + c 3 )a 2 ] , 

where b := a z , ip := <p z , and 7 := 1 + c\ . Solutions of the ODEs (||) correspond to coherent 
structures of the CGLE. 

The simplest relevant solutions of these ODEs are the fixed points given by (a, b, ip) = 
(a/1 — q 2 , 0, q); these correspond to plane wave solutions of the CGLE where A{x, t) = yl — q 2 exp i(qx + ujt) 
and uj = C3 — q 2 (ci + C3). An example of more complex solutions of the ODEs (||) are heteroclinic 
orbits which correspond to coherent structures that asymptotically connect different states. Ex- 
amples of such structures are fronts that connect nonlinear plane waves to the homogeneous 



state A = [24] and Nozaki-Bekki holes that connect plane waves of different wavenumber q 
|24j,|9|. 

Here we present an extensive study of the structures that are associated with the limit cycles 
of the ODEs (||) |3(|. These limit cycles correspond to spatially periodic solutions of the CGLE 
that we have already referred to as MAWs (Fig. ||). For appropriate choices of c\ and C3, the 
period P of these MAWs can be made arbitrarily large, and in this limit the limit cycles approach 
a homoclinic orbit connecting the stable and unstable manifold of one of the plane wave fixed 
points (Fig. ||a). Some of these infinite period MAWs have also been referred to as "homoclinic" 
holes, and have been studied extensively recently [25[,31|; they are qualitatively different from 
the well-known Nozaki-Bekki holes 



Even if the coefficients c\ and C3 are fixed, MAWs are not uniquely determined. Counting ar- 
guments, similar to those developed in [p4]| , yield that in general we may expect a two-parameter 
family of solutions. Let us first perform the counting for the homoclinic orbits. As shown in 
|2q| , these orbits connect the one-dimensional unstable manifold of a fixed point with its two- 
dimensional stable manifold. In general, one needs to satisfy one condition to make such a 
connection, in other words, such a homoclinic orbit is of codimension one. Since the coherent 
structure Ansatz (^) has two freely adjustable parameters {uj and v), we therefore expect a one 
parameter family of homoclinic orbits. 

The situation for the limit cycles of the ODEs is even simpler. Limit cycles are of codimension 
zero in parameter space, and so we expect a two parameter family of limit cycles. In other words, 
if we have found a limit cycle for certain values of v and uj, then we expect this limit cycle to 
persist for nearby values of the parameters v and uj. 

Obviously, we can parameterize this family of limit cycle coherent structures by v and uj, but 
this is not very insightful. Instead we will use the following two quantities that are more directly 
accessible in studies of the CGLE: the spatial period P of the MAWs, and their average phase 
gradient v := (Jq dx<p x )/ P. Note that for homoclinic holes, P simply goes to infinity; thus 
homoclinic orbits and limit cycles are members of a single family. 
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The multiplicity of the MAWs can also be obtained by considering the instability of the plane 
wave solutions from which the MAWs emerge [^] (see section |2.2| below). The plane waves 
form a one-parameter (q) family and undergo the well-known Eckhaus instability when the 
coefficients ci, C3 are increased beyond certain critical values which depend on q. In the unstable 
regime, a plane wave with wavenumber q is unstable to a whole band of perturbations with 
wavenumbers k € [0, k max (q)] j|. For finite systems of size L, this instability thus only appears 
when L > L m i n = 2ir/k max . Therefore for each q there is a unique one-parameter (L) family 
of perturbations that can render the plane wave unstable and at each of the corresponding 
bifurcations a new MAW solution emerges. Hence also by this line of reasoning MAWs form a 
two-parameter family. 

2.2. Bifurcation Scenario for MAWs 
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Figure 4. Bifurcation diagrams for fixed C3 = 2 and P = 30, showing Hopf (filled square), 
drift pitchfork (open diamond) and saddle-node (triangle) bifurcations. The dot-dashed line 
represents the homogeneously oscillating solution of the CGLE, while lower and upper branch 
MAWs are represented by full and dashed curves respectively, (a) Overview of the maximum 
phase gradient of the MAWs as function of c\, (b) close-up, (c) the minimum of \A\, and (d) the 
velocity v. For details see text. 



The general counting arguments given in the previous section do not provide information on 
the range of existence of MAWs as a function of the coefficients c\ and C3 and the parameters 
v and P. Here we will focus our analysis on the v = case since this is most relevant for the 
transition to defect chaos [33]; the v ^ case will be treated elsewhere p4| . 

All bifurcation computations have been performed with the aid of the software package 
AUT094 [35|. AUT094 can trace MAW solutions through parameter space, and when it detects 
bifurcations it can follow the newly emerging branches. AUT094 discretizes the ODEs (§) on 
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Figure 5. (a) Phase gradient and (b) amplitude profiles of a lower branch (full curve) and upper 
branch (dashed curve) MAW, obtained for c\ =0.6, C3 = 2, P = 30. 



a periodic domain of length L, and L will play the role of the period P of the MAWs. Control 
of the average phase gradient v = i>q is implemented via the integral constraint J tpdz = Luq. 
Since periodic boundary conditions result in translational invariance, we introduce an additional 
"pinning" condition a z (0) = in order to obtain unique solutions. 

Under these conditions, the continuation procedure works as follows. First of all, v and P are 
set to fixed values, and throughout this paper we will set v = 0. Starting from a known solution 
such as a plane wave or a coherent structure obtained by other means, AUT094 is set up to 
trace the MAWs along trajectories in ci, C3 space, while calculating the parameters uj and v of 
these MAWs. 

The results of our bifurcation analysis are summarized in Fig. [|. When c\ or C3 is increased, the 
uniformly oscillating state of the CGLE (A(x, t) = e lC3t ) becomes unstable via a Hopf bifurcation, 
from which stationary MAWs emerge (section |2.2.1|) . These stationary, left-right symmetric 
solutions undergo a drift pitchfork bifurcation, which leads to left and right traveling MAWs 
(section [2,2,2 , see also Fig. ||b); as discussed later, these are the solutions relevant for the 



dynamics in the phase chaotic regime. Following these branches of traveling MAWs, we encounter 



a saddle-node bifurcation where an "upper" and "lower" branch of MAWs merge (section |2.2.3 



see also Fig. ||); this bifurcation limits the range of existence of MAWs and is closely related to 
the formation of defects. The upper branch MAWs can be continued back to negative values of 
ci, where they terminate in a solution consisting of a periodic array of shocks and stationary 
Nozaki-Bekki holes [29). Upper branch MAWs with P — > 00 have been studied under the name 



homoclinic holes 25, 3j 



It should be noted that, without loss of generality, we focus here on solutions with v > 0, for 
which the main peak of the phase gradient profile is positive (see Fig.|5|). Solutions with v < 
can be obtained from right moving MAWs by applying the mapping x — > —x, z — > —z, v — > 

-v,a z -> -a z ,(p z -> -<f> z . 

2.2.1. Benjamin-Feir instability - Hopf bifurcation 

Since the average phase gradient v is conserved across bifurcations, we start the continuation 
procedure from the uniformly oscillating solution A(x, t) = e tcst that has v = 0. On an infinite 
domain this uniformly oscillating solution becomes unstable via the so-called Benjamin-Feir 
instability when C1C3 > 1 ||. In a finite domain of size L, the onset of this instability is shifted 
to higher values of the product C1C3 [36|; this finite size effect is relevant for our studies since 
the spatial period L = P is fixed in the continuation procedure. 
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Figure 6. (a) Location of Hopf (dashed curve) and drift pitchfork bifurcation (dotted curve) 
in cx,P space for 03 = 2.0. (b) The shaded area reported in the ci,cs coefficient space indicates 
where the drift pitchfork bifurcation does occur. The thick full curve in (a,b) indicates the 
Benjamin- Feir-Newell instability for infinite domains, (c) Example of a bifurcation diagram for 
large values of the coefficients c\ = 10, C3 = 5 where the drift pitchfork bifurcation does not 
occur. For increasing P the MAW solutions approach regular arrays of stationary pulses; an 
example of such a pulse is shown in (d) for P = 30. 



In the ODEs (Q), the fixed point (a, 6, ip) = (1,0,0) corresponds to the homogeneously oscil- 
lating solution. For given values of the period P, this fixed point undergoes a Hopf bifurcation 
at values of c\ and C3 where in the CGLE (|l]) the mode with wavenumber 2n/P becomes unsta- 
ble |56|. This Hopf bifurcation was analytically shown to be supercritical for sufficiently small v 
and large P in earlier studies [p!^j32]| ; our numerical results are consistent with this. For finite P, 
the solution bifurcating from the fixed point is a limit cycle which approaches a homoclinic orbit 
in the limit P — > 00. The solutions of the CGLE that correspond to these orbits are stationary, 
reflection symmetric MAWs; an example of these is shown in Fig. ^a. 

2.2.2. Drift pitchfork bifurcation 

When the CGLE coefficients c\ and/or C3 are increased further, the stationary MAW under- 
goes a drift pitchfork bifurcation |i7j from which two new branches of asymmetric (v 7^ 0) MAWs 
emerge (see Fig. |I|b); one of these moves to the left, one to the right. The locations of both 
the Hopf and the drift pitchfork bifurcation approach the Benjamin- Feir-Newell curve for large 
P (Fig. ||a), while for smaller P the drift pitchfork occurs for increasingly larger coefficients 
c\ and C3. However, only when these coefficient lie in the range shown as the shaded area in 
Fig. |6]b, the pitchfork bifurcation can occur. Otherwise, only stationary MAWs are found. For 
increasing c\ and C3 these MAWs become pulse-like and finally approach the solitonic solutions 
of the nonlinear Schrodinger equation |24j (Fig. oc,d). 

For the case v 7^ [|34}j, the initial plane wave already breaks the reflection symmetry, the 
initial MAW has nonzero velocity and the drift pitchfork bifurcation is replaced by its typical 
unfolding p8|. 
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2.2.3. Saddle-node bifurcation 

Along the branch of right traveling MAWs that we described above, the maximum of the 
phase gradient grows with increasing c\ and C3 until a saddle-node (SN) bifurcation is reached, 
where these MAWs merge with another branch of MAW-like solutions. To distinguish these 
branches we refer to them as the "lower" and the "upper" branch; for examples see Figs. ^J8|. 
The lower branch MAWs are the key to understand more of the phenomenology of phase chaos. 
The upper branch MAWs can, similarly to the lower branch MAWs, be parameterized by v and 
P, but for the same parameters, they present more pronounced modulations (see Fig. ^). 

The most important aspect of the saddle-node bifurcation is that it limits the range of existence 
of MAWs, since we will show that this limit is responsible for the transition from phase to defect 
chaos. Fixing i/ = 0, the locations of these bifurcations form a two-dimensional manifold in the 
three dimensional space spanned by c±, C3 and P. In Fig. [7|a the saddle-node curves are shown 
in the 01,03 coefficient plane for a number of fixed periods P; for larger P, the values of ci,C3 
where the bifurcation takes place decrease. In Fig. ^b the saddle-node curves for a number of 
fixed values of C3 are shown in the P, c\ plane; for larger C3 (ci), the saddle- node occurs for 
smaller values of P and c\ (03) Once the coefficients c\ and C3 are fixed, we define Psn 

as the period for which the saddle-node bifurcation occurs. Note that there is also a range of 
coefficients c\ and C3 (between the P — > 00 and c\cz = 1 curve where the saddle- node bifurcation 
does not occur. 



2.3. Evolution of perturbed MAWs 

In this section we will show that many basic aspects of the phenomenology of the CGLE can be 
understood from a typical bifurcation diagram of MAWs such as shown in Fig. ||. We have chosen 
fixed coefficients c\ = 0.65 and C3 = 2 and varied the spatial period P of MAWs that exist at 
these coefficients. Three families of solutions are represented: the homogeneous oscillation, the 
lower branch (LB) and the upper branch (UB) MAWs. The shaded area schematically indicates 
the near-MAW structures observed in phase chaotic states such as shown in Fig. |l| (f,i,j). The 
arrows in Fig. ^ represent the dynamical evolution of perturbed MAWs, and their direction can 
be obtained by performing a linear stability analysis. 

Linear stability - As discussed in section the homogeneous solution is stable against 
short wavelength perturbations (arrow 1), and turns unstable via the Hopf bifurcation that also 
generates the lower branch MAWs (arrows 2). As discussed in 1 25, 31 1 , upper branch MAWs 
have at least one unstable eigenvalue, and the dynamical evolution of perturbations is directed 
away from upper branch MAWs (arrows 3,4). 



11 




interpeak spacing p P SN =52.6 

Figure 8. Illustration of the relation between MAW structures and phase chaos for c\ = 0.65, C3 = 
2. Both for MAWs and for an extended profile obtained from a phase chaotic state one can 
extract the values of subsequent phase gradient peaks (vertical axis) and their inter-peak spacing 
(horizontal axis). The curves show the bifurcation diagram for the lower branch (LB) MAWs 
and upper branch (UB) MAWs while the shaded area indicates the typical values for near-MAW 
structures that occur in phase chaos. Full (dashed) curves denote stable (unstable) solutions for 
system size L = p. Arrows show the typical evolutions of near- MAWs. The coefficients here are 
equal to those in Fig. |l| (f,i,j) and for this case we have found that phase chaos is only a long 
lived transient: the shaded area reaches Psn just before defects appear. 
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Figure 9. The evolution of an unstable homogeneous state towards lower branch MAW dynamics, 
for c\ = 3 and C3 = 0.6. The coefficients c\ and C3 are chosen such that no saddle-node bifurcation 
occurs for any P. (a) Evolution towards a stable stationary lower branch MAW for system size 
L = 25 and (b) towards a stable drifting lower branch MAW for system size L = 30. Note that 
for the coefficients chosen, the drift pitchfork bifurcation occurs at P = 27.7. (c) Evolution 
towards phase chaos for system size L = 100. Incoherent evolution of structures characterized 
by local concentrations of phase gradients can be clearly observed. We think of these structures 
as "near" MAWs. 
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Figure 10. Evolution of perturbations of the upper branch MAWs. The coefficients c\ and C3 are 
chosen such that no saddle-node bifurcation occurs for any value of P. (a) A slowing down and 
spreading of the phase gradient characterizes the decay to a lower branch MAW for c\ = 0.55 
and C3 = 2. (b) For the same coefficients, another perturbation leads to an increase in velocity 
and divergence of the phase gradient. A defect occurs, from which hole-defect dynamics spreads 
for these coefficients, (c) For c\ = 3 and C3 = 0.6 a perturbed upper branch MAW leads to a 
defect, but defects do not percolate through the system. 



The linear stability of lower branch MAWs will be discussed in more detail in section ^J. It 
turns out that perturbations of lower branch MAWs can evolve in many ways, but in almost 
all cases the ensuing dynamics remains close to the lower MAW branch (shaded area in Fig. |8|). 
The only exception we have found to this rule is when a MAW is pushed beyond the saddle-node 
bifurcation (arrow 5). 

Nonlinear evolution - Here we want to go beyond the linear analysis and study the nonlinear 
evolution of MAWs along the arrows of Fig. ||. The examples (at different choices of the coef- 
ficients) of the dynamics shown below are not exhaustive, but should serve to illustrate typical 
behavior which appears to be very robust. 

arrow 2 - When the uniform oscillation becomes linearly unstable perturbations grow. To 
the left of the saddle-node, perturbations evolve to dynamics dominated by lower branch MAWs 
(Fig. |9|). For small system sizes, stable MAWs may occur (Fig. ^a,b), while for larger systems 
periodic sequences of MAWs are unstable with respect to the so-called interaction or splitting 
instabilities (HU^I that will be discussed in section 3.4. Hence a perturbed unstable homogeneous 
state typically does not converge to a train of coherent MAWs, but instead evolves to phase chaos 
(Fig. ||c) . In the context of the bifurcation diagram, note that the disordered structures observed 
in the phase chaotic evolution are quite similar to lower branch MAWs. The shaded area in Fig. |8| 
represents this "near-MAW" behavior. 

arrows 3,4 - Upper branch MAWs are always unstable due to the positive eigenvalue asso- 
ciated with the saddle-node bifurcation. The resulting incoherent dynamics has been studied 
quite extensively in the context of hole-defect dynamics [25,31]. (i) When a perturbation has 
pushed an upper branch MAW towards the "lower" part of the bifurcation diagram, the struc- 
ture decays towards lower branch MAWs (arrow 3). An example of a space time plot for the 
decay towards a lower branch MAW is shown in Fig. [l0|a. (ii) When the perturbation pushes 
the MAW towards the "upper" side of the diagram, the phase gradient peak that characterizes 
MAWs grows without bound, and at the same time the minimum of \A\ approaches zero: a 
defect is formed (arrow 4). The dynamics after such a defect has formed depends on the values 
of the coefficients c\ and C3. Two different examples are shown in Fig. [l0[b,c. For more details 
see section || 

arrow 5 - So far we have encountered two scenarios: if the phase gradient peak of a structure 
is "larger" than that of an upper branch MAW, then it will grow out to form defects. If it is 
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"smaller", it will decay back in the direction of the lower branch MAWs. The latter process 
frequently occurs in phase chaos, preventing the formation of defects, while the former process 
needs to be initiated by appropriate initial conditions. However, when the upper and lower 
branches approach each other and disappear in a saddle- node bifurcation, there are no structures 
left to prevent arbitrary small perturbations to grow out to defects. This dynamical process, 
which is represented by arrow 5 in Fig. ||, is the core of our argument: defect formation takes 
place beyond the saddle-node bifurcation. 

2.4. Breakdown of phase description 

An alternative approach to describe the creation of defects from phase chaotic states is via 



blow-ups in so-called phase-equations [2C]. Phase equations are based on the observation that 
close to the onset of phase chaos (near the Benjamin- Feir-Newell curve) the amplitude is "slaved" 
to the phase dynamics. In this situation a phase equation can be obtained by a gradient expan- 
sion Q. The expansion including all parity-symmetric terms up to fourth order |2C] reads 

+ ^ 2) di&± + ( 3) /W\ 2 + (4) m\ 2 &± 

4 dx dx 3 4 I dx 2 I 4 V dx ) dx 2 



where 0^ = 1 - cic 3 , ^ 2) =-(<* + c 3 ), djp = -c\ (1 + c|)/2, = -2ci(l + c|), df = -ci(l + 
c|), Q^=— 2(1+03). ^he l° wes t order description of phase chaos is obtained when the parameters 
nf\nf } and are set equal to zero; the resulting equation is known as the Kuramoto- 
Sivashinsky equation [§| . 

The phase equations with higher order terms included have been studied via direct integration 
by Sakaguchi pCp . For the full Eq. (0), Sakaguchi observed finite time divergences of the phase 
gradient for coefficients close to the transition from phase to defect chaos in the CGLE. He 
attributed such divergences to the occurrence of defects in the CGLE. No blow-up of the phase 
gradient is observed for Eq. (0) without the last term, or for the simple Kuramoto-Sivashinsky 



equation. Recently, Abel et al. [41] quantified the increasing discrepancies between the phase 
equations of different orders and the full dynamics in the CGLE with increasing distance from 
the Benjamin- Feir-Newell curve and identified the relative importance of the various terms in 
Eq. (1). 

Since the essential ingredient of our theory is the occurrence of a saddle- node bifurcation, we 
have investigated the bifurcation scenario for various truncations of the phase equations (||). In 
the context of phase dynamics, our Ansatz (g) becomes of the form 

4>(x, t) = 4>{x — vt) + (to — Cs)t . (5) 

We have studied MAW-like structures occurring in the phase equations by employing the same 
methodology as for the CGLE; the average phase gradient value v is fixed to and P parame- 



terizes the spatial period of the MAW. In Fig. 11 we compare bifurcation diagrams and MAW 



profiles for different expansions at the parameters c\ = 3.5, P = 50. 

For all phase equations considered here the coherent structures are again born in a Hopf and 
undergo a drift pitchfork bifurcation, beyond which the maximal phase gradients increase. This 
leads to increasing discrepancies between different approximations. In particular, the coherent 
structures for Eq. (|j) exhibit saddle-node bifurcations at parameter values not far from those 
for corresponding MAWs in the CGLE; nevertheless the MAWs of Eq. (0) deviate substantially 
from the CGLE MAWs for the upper branch of MAWs. The Kuramoto-Sivashinsky equation, 



14 




0.4 



0.2 



-0.2 



1 1 1 

: A 


1 1 1 1 


1,1,1,1, 



10 20 x 30 40 50 



Figure 11. Comparison of different phase expansions : Eq.(|) (dotted), Eq.(§) without the last 
term (dot-dashed), Kuramoto-Sivashinsky Eq. (dashed) and CGLE (full curve). Parameters 
are c\ = 3.5, v = 0, P = 50. (a) Bifurcation diagrams and (b) spatial profiles of lower branch 
coherent structures at C3 = 0.7. 



and Eq. (Q) without the last term, do not exhibit a saddle-node bifurcation. Since these latter 
two models do not experience blow-up, we can safely conclude that these observations confirm 
our picture, and that the saddle-node bifurcations of coherent structures play the same crucial 
role in both the full CGLE and its phase equations. 



3. Large scale chaos 

In this section we will study the dynamical evolution of the CGLE near the transitions from 
phase to defect chaos. The transition between these two states can either be hysteretic or 
continuous: in the former case, the transition is referred to as L3, in the latter as L\ | |42|] . 

How are defects generated from phase chaos? Let us start to consider a small system in which 
a stable lower branch MAWs has been created. When we fix the coefficients c\ and C3 and 
steadily increase the size of the system, and hence the period P of the MAW, we find that as 
soon as we push P beyond Psn, the MAW structure blows up to form defects. An example of 
this is shown in Fig. [l^a. In a similar fashion, defects are created when the system size L is 
fixed, and either c\ or C3 are increased until Psn<L (Fig. |]c,d). 

How is this related to phase chaos? As shown in Fig. |T2|b, typical phase chaotic states show 
much more incoherent dynamics, containing many MAW like structures but of much smaller 
period. Our central conjecture is therefore that the transitions from phase to defect chaos are 
triggered by the occurrence of near-MAW structures in a phase chaotic state with v = |33| and 
periods larger than Psn, the spatial period of the critical nucleus for defect creation. 

To test this conjecture, we have numerically investigated the distribution of inter-peak spacings 
p of the phase gradient profile (see Fig. [j]e,f ). In section 34. we discuss the definition of p and the 



details of our numerical analysis. In particular, we have examined in the c\, C3 plane 17 different 



"cuts" across the L\ and L3 transition lines. In section 3.2 the results of our numerics along a cut 



through the L\ transition line are presented, while section |3.3| is devoted to the L3 transition. 



We will show that the presence of inter-peak spacings p larger than Psn accurately predicts 



the transition from phase to defect chaos (Fig. |2|). In the last section 3.4 we will show that a 
reasonable, parameter-free estimate of the numerically observed transitions can be obtained via 
a linear stability analysis of the MAWs. 
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Figure 12. Defect formation at C3 = 2, c\ = 0.7. (a) Defect formation. As initial condition we 
took a lower branch MAW with P < Psn = 26.8 which we embedded in a background of zero 
wavenumber. The system size L here is equal to 28.8, which is larger than Psn and a defect is 
formed; for L < P$n this defect formation does not take place, (b) Random initial conditions in 
general evolve to MAW like structures with P < Psn which do not lead to defects; the "critical" 
nucleus that leads to defect formation has a rather small basin of attraction here. 



3.1. Identification of MAWs in the phase-chaotic regime 

To verify our main conjecture, we have to characterize the MAW structures occurring in the 
phase-chaotic regime. In general this is a complicated task, since the phase gradient profile 
of a typical phase chaotic state (see Fig. |l]e,f and [b|) consists of many peaks of different size, 
spacing and shape; a priori it is unclear how to compare these to MAW profiles. However, a 
close inspection of the defect forming process reveals that while closely spaced phase gradient 
peaks evolve in a quite erratic way, well spaced peaks appear to have a more regular dynamics 
and frequently their overall shape resembles that of MAWs (see Fig. 13). These large period 
near-MAWs modify their shape quite slowly with respect to the other structures present in the 
chaotic field, and propagate over a disordered background. Therefore we study the distribution 
of inter-peak distances p, keeping in mind that the tail of this distribution is relevant for defect 
generation. 

The phase gradient profile of a coherent MAW (see Figs. [l|a and Fig. |E|a) shows a secondary 
maximum. To obtain the correct period P of a near-MAW, such small extrema should be 
neglected when the inter-peak spacing p is measured. We introduce a cutoff for the size of the 
phase gradient peak equal to the size of the secondary extremum of the MAW with the largest 
P. As an additional result of this cutoff, small fluctuations are not considered as MAW peaks. 
It should be noted that the tail of the distribution of p is rather insensitive to the precise value 
of this cutoff. 

In order to estimate the probability density D(p), for every time interval r = 0.5, the inter- 
peak periods p of the spatial profile of the phase gradient are determined. In addition, for every 
snapshot the largest value p max of the inter-peak spacing p is stored separately, and this leads 
to the distribution D(p max ). From the spatial profile of \A\ the distribution and the 

minimal amplitude value |^4| m i n can be derived. This latter quantity is used to detect defects: 
when |j4| m j n falls below a value of 0.1, we take this as an indication of a defect. 

Extensive simulations have been made possible thanks to an innovative time-splitting code 
which ensures precision and stability comparable with pseudo-spectral codes, but is noticeably 
faster [23]. The spatial resolution Ax has been set to 0.5 and the integration time step to 
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Figure 13. Local phase gradient of the chaotic field just before defect formation for c\ = 0.65 
and C3 = 2. Panel (a) is a snapshot of the field at a time t = 120 before the occurrence of the 
defect, (b) , (c) and (d) are successive snapshots taken at time intervals St = 30. In (a) also the 
shape of the MAW at the saddle-node is superimposed (thick dashed line) on the profile. 



0.05. Simulations have been carried out for integration times ranging from t 



10 & to 



t = 3 x 10 7 and for a typical system size L = 512; occasionally, runs have been performed with 
L = 100, 200 and 5000. Typically, our runs start from random initial conditions of the type 
A k (t = 0) = \A\ k (t = 0) • ^^=0) ( w here A k (t) = A(kAx,t) and <f> k (t) = <f>(kAx,t)) with 



,(0) 



\A\ k (0) = l + r k 
fc _!(0)*0.8 + % 



(6) 
(7) 



where r k and q k are random numbers uniformly distributed in [—0.05, +0.05] and </>i(0) = 0.0005. 
This initial condition (|7|) leads to a smooth phase and the formation of defects due to initial 
discontinuities is avoided. 

In sections ^2 and we will consider in detail two particular cuts in the (01,03) coefficient 



space, one across the L\ and one across the L3 curve. In particular, we will analyze the behavior 
of the probability densities Z)(|A|), D(jp) and D(p max ) for both transitions. 

3.2. L\ transition 

In this section we concentrate on the L\ transition that is observed when the value of c\ is 
fixed at 3.0 and C3 is varied. 

Transition to defect chaos - Starting from random initial conditions we have integrated the 
dynamics of the CGLE for long durations. For a fixed system size L we observe that, as a 
function of the total integration time, the value of c| for which defects are formed appears to 
decrease. Similar behavior occurs when the system size L is increased for fixed integration times. 
For example, for an integration time of 3 x 10 and c\ = 3 we find for system size 100, 200 and 
512 critical values 0.82,0.81 and 0.79, respectively. For a size L = 5000 and integration times 
3 x 10 6 a critical value of 0.79 is also found. 
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Note that even the lowest value of C3 for the numerically measured transition obtained here 
is far above the lower bound C3 = 0.704 which is the value of C3 where the size of the critical 
nucleus for defect formation diverges (Psn —* 00). Below, we will give an estimate of the critical 
value 63 for which the defect density should vanish in the thermodynamic limit by extrapolating 
finite time and finite size data. 

Distribution ofp - Let us now consider the distribution of p's for various coefficients C3 near the 
L\ transition. It is clear from the data reported in Fig. [14| that the shape of these distributions 
is quite insensitive to the presence or absence of defects. This can be partly explained by 
the fact that just above the L\ transition defects arise in the system as rare isolated events 



occurring during the spatio-temporal evolution, as shown in Fig. 10:. This is fully consistent 



with earlier observations that the L\ transition is continuous 17,0| 21 1. We focus on the tail of 
the probability density D(p), since this gives information on the probability to observe defects. 
Our numerical results suggest an exponential decay, i.e., D(p) oc exp(— a ■ p) with a = 0.6 for 
sufficiently large p. 

Similarly to the apparent transition value c|, the values associated to extremal events |A| m j n 
and p m ax depend on integration times and system sizes. By assuming that D{p) remains finite 
(but likely exponentially small) for large p, we can expect that for long enough times, rare events 
associated with large values p will occur, and hence, defects can form after possibly very long 
transients. 

Crossover behavior - A good order parameter to identify the occurrence of the transition 
starting from the defect chaos phase near the L\ transition is the defect density <5d which 
measures the number of defects occurring per space and time unity. In the defect chaos regime 
80 > 0, while it vanishes at the Li-transition. Now we can relate this order parameter to the 
tail of the distribution of p. Our conjecture states that defects should arise when p > Psn, 
therefore the defect density 5n should be related to the probability to have structures of period 
p > P S n, i-e., 

POO 

5 d oc / dpD(p) oc e ~ aPsN ; (8) 



where D(p) oc exp(— a ■ p) has been used. If we now assume that the distribution D(p) does 
not vary significantly across the transition (as is evident from Fig. |l4|), then the change in the 
probability to have p > Psn is dominated by the changes in Psn with C3. A reasonable fit of 
our bifurcation data for Psn (see Fig. f^) in the interval 30 < Psn < 300 is 

Psn ~ — ^ , (9) 
c 3 - cf 

where (3 ~ 4.38. Combining this result with the Ansatz (||), we immediately obtain the following 
expression for the defect density: 

5 d oc e -^/( c 3- c 3°) . (10) 

A similar expression was proposed in ]l7| , f2l"| for the defect density near the L\ transition. 

In order to verify if the expression ([ji]) is reasonable also for our choice of the parameters, we 
have estimated the probability |23| 

w(\A\)= r d\A\D(\A\) , (11) 
J 

to observe an amplitude less than \A\. This quantity gives a more precise characterization of the 
Li-transition than 5d, because it measures not only the extreme events corresponding to true 
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Figure 14. Probability densities (a) D(p) and (b) D(p max ) for c\ = 3 and various values of C3 
on a lin-log scale. The curves refer to C3 below C3 = 0.79 (namely to C3 = 0.77 and 0.78), as well 
as to values corresponding to the defect chaotic regime: C3 = 0.79, 0.80 and 0.81. The system 
size was L = 512 and the integration times where t = 5 x 10 5 for C3 = 0.81, t = 5 x 10 6 for 
C3 = 0.80 and t = 25 x 10 6 for all other values. 
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Fi gure 15. Probability densities (a) D{p) and (b) D{p max ) for C3 = 2 and various values of c\ 
reported in a lin-log scale. The data are for a system size L = 512 and for integration times 
ranging from t = 5 x 10 5 for c\ = 0.63,0.64,0.66 and 0.65 (PC)* to t = 2.5 x 10 6 for ci = 0.65 
(DC). The labels DC and PC indicate that we are in presence or absence of defects, respectively. 
The label (PC)* refers to the regime before defect formation at c\ = 0.65. 



defects, but also the tendency of the system to generate structures characterized by small |^4| m j n - 
We estimated the quantity ([ll]) for several \A\ values and for various C3 parameter values in the 
defect chaos regime. Reporting ln[u)(|^4|)] as a function of 1/(03 — £3) a reasonable linear scaling 
is observed in the range 0.795 < C3 < 0.85, for 0.1 < \A\ < 0.5, with the choice C3 = 0.72. The 
value C3 where the defect density should asymptotically vanish is much smaller than c\ obtained 
via direct numerical simulations but still bigger than C3 = 0.704 where Psn — ¥ 00. 

We can now easily estimate the integration time needed to observe a tiny shift of the apparent 
value C3 towards the corresponding asymptotic value C3 ~ 0.704. Limiting our analysis to system 
size L = 512, a typical time-scale to observe a defect at C3 = 0.79 is t ~ 3 x 10 7 . At this value 
of C3, P5at = 46.5, while for C3 = 0.739, Pstv = 105 • Invoking the exponential decay of D(p), one 
immediately finds that the time scale to observe a defect at C3 = 0.739 is of order 10 17 , which is 
completely outside the reach of present day computers. 
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Figure 16. Results of the linear stability analysis : (a) leading part of the eigenvalue spectrum 
(continuous spectrum denoted by the full curve, Goldstone modes by a filled square, saddle- 
node by filled triangle, interaction by open square and splitting modes in L = P by dots, 
respectively), (b) splitting eigenmodes (dot-dashed and dotted, \ S pUt = 0.018 ± 0.28i) of the 
phase in L = P compared with spatial MAW profile of the phase gradient (full curve). (c,d) 
Space time plots showing the splitting of a MAW initially perturbed by small noise. Parameters 
are c\ = 3, C3 = 0.72, P = 43 near L\ for (a-c) and c\ = 0.65, C3 = 2, P = 35 near L3 for (d). 



3.3. L3 transition 

In order to characterize the L3 transition from phase to defect chaos in more detail C3 = 2 



has been fixed, while the coefficient c\ is varied. The L3 transition is hysteretic |17,l9fl: to the 
left of L3 one may have phase or defect chaos depending on the initial conditions. Beyond the 
L3 phase chaos breaks down and defects occur spontaneously for any initial condition. In order 
to study the dynamics across this transition we therefore initialized the simulations with initial 
conditions ©,(0) or used relaxed phase chaos configurations corresponding to values of c\ far 
below the L3 line. 



The probability densities D{p) and D{p max ) are shown in Fig. 15. For c\ < c* = 0.65 all 
distributions collapse on a unique curve, but as soon as defects arise the distributions change 
substantially. Whenever a defect is generated, hole-defect dynamics takes place (see Fig. |i~0|b). 
As a result phase chaos is replaced by defect chaos. The noticeable modification of the distribu- 
tions thus reflects the fact that the L3 transition is discontinuous. Also the probability density 
for \A\ changes abruptly across the L3 transition. 

3.4. Mechanism for the selection of p 

When approaching the transition to defect chaos from the Benjamin- Feir-Newell curve, three 
parameter regions, corresponding to different dynamical regimes, can be distinguished (Fig. |i~8|) . 
The first encountered region corresponds to infinite values of Psn '■ here we expect no defects 
to occur, irrespectively of system size and integration time. The phase chaos is the asymptotic 
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Figure 17. Results of the linear stability analysis : (a) leading part of the eigenvalue spectrum 
(continuous spectrum denoted by the full curve, Goldstone modes by a filled square, saddle-node 
by filled triangle and interaction by open square at Xi n t = +0.0048, respectively), (b) Goldstone 
mode (dot-dashed) and interaction eigenmode (dotted curve) for the phase compared with spatial 
MAW profile of the phase gradient (full curve) in L = 2P. (c,d) Space time plots showing the 
attraction of two periods of the same MAW initially perturbed by the interaction eigenmode. 
Parameters are c\ = 0.7, C3 = 2, P = 25 near L3 for (a-c) and c\ = 3, C3 = 0.85, P = 21 near L\ 
for (d). 



regime in this first region. Then, when c\ and/or C3 are increased, a crossover regime is reached 
where extreme events (large inter-peak spacings) may lead to defect formation. Here phase 
chaos can persist as a long lived transient, but eventually we expect it to break down. Then, 
when ci and/or C3 are even further increased, we experience a dramatic drop in transient times, 
and defect chaos sets in quite rapidly. We understand this drop to occur when typical values of 
p (and not rare extreme events) become larger than the corresponding P$n values. 

An approximate prediction for the location of the apparent phase to defect chaos transition 
(numerically obtained from the defect density) can be achieved in terms of a simple linear 
stability analysis of the MAWs (Figs. |l6| and |l7|). A key element in our framework is the 
"typical large value" of p as a function of coefficients c\ and C3; below we will identify two linear 
instabilities that act to either increase or decrease p, and their balance sets a scale for typical p 
that will predict the location of the transition from phase to defect chaos rather well. 

Due to translational and phase symmetries both MAW branches have neutral modes, i.e., 
Goldstone modes. The eigenvalue associated with the saddle-node bifurcation is positive for 
MAWs of the upper branch and negative for the lower branch. In what follows the lower branch 
MAWs are considered exclusively. 

Splitting - The spatial structure of a MAW of large period consists, roughly, of a homogeneous 
plane wave part and a local peak part. For the parameter regime we consider here, fully extended 
plane waves are linearly unstable, and so we may expect that the MAW spectrum will be 
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dominated by this instability for sufficiently large values of P. Our linear stability analysis 
indeed shows that for appropriate parameters {L = P) and small enough P, all eigenvalues 
Aj < 0, but when we increase P, MAWs become linearly unstable (X S plit > 0> Fig. 16). The 
shape of the unstable eigenmodes (Fig. |l6|b) suggests that this instability leads to the growth 
of a new peak in the homogeneous part of the MAW, and this is indeed the behavior observed 
in numerical simulations of the perturbed MAW (Fig. |i~6|c,d). As a result two (or more) short 
MAWs with smaller P will appear. We interpret this process as the splitting of a MAW in 
two or more smaller MAWs and we call the eigenmodes associated to such instability "splitting 
modes" . 

Clearly, this instability tends to reduce the peak-to-peak distances p and prevents MAWs 
to cross the SN boundary; in the phase chaotic regime this instability tends to inhibit defect 
generation. 

Interaction - By using a Bloch Ansatz [ 43 1 , we extended the stability analysis to systems 
with n identical pulses (L = nP). For n > 1, an additional instability may appear [Q] (see 
Fig. |l7| ). Eigenvalues \i n t > are found mainly for small P (typically P < 30). The shape of 
the eigenmodes, i.e., an alternating sequence of positive and negative translational Goldstone 
modes (Fig. |i~7|b), suggests that the instability is due to the interaction between adjacent MAWs. 
This interaction shifts adjacent peaks into opposite directions, thereby creating occasional larger 
values of p (Fig. |i"7|c,d). In phase chaos this process leads to an increase of the spacing p between 
some peaks, thus enhancing the generation of defects. 

Competition of Instabilities - Both the splitting and interaction mechanisms are similar to 
instabilities observed in the Kuramoto-Sivashinsky equation . We believe that phase chaos 
is governed by the competition of these two mechanisms that tend to increase or decrease the 
inter-peak spacings p. Almost independent of the coefficients the splitting instability dominates 
for MAWs with P > 30. This can explain why large inter-peak spacings p > 30 become rare as 



reported in Figs. 14,15 



We suggest a connection between the interchanging dominance of these two different insta- 
bilities and the sudden change of 5d (near L\) or the transient times before defect occurrence 
(near L3). We calculated the linear stability spectra for a variety of coefficients and periods P 
close to Psn- From these we obtain a curve in coefficient space (Fig. |jl]) where the real parts of 
interaction and splitting eigenvalues are equal. For larger c\ or C3, interaction becomes stronger, 
and we expect larger p 7 s and defect formation, while for smaller c\ and C3, splitting dominates, 
p's are decreased and defect formation becomes rare. 

As shown in Fig. 18, the curve where the two instabilities are equally strong near the saddle- 
node bifurcation gives a rather good estimate of where the apparent transition from phase to 
defect chaos occurs. Notice that in this "balance of instabilities" picture, there is no tunable 
parameter: once we have calculated Psn and the instabilities of the MAWs for a range of 
coefficients, a precise prediction for the "transition" from phase to defect chaos can be given. 



4. Discussion and Final Remarks 

In this section we report some open questions related to defect formation, together with some 
final remarks and a brief outlook. 

Further Refinements - In order to accurately test our results, we have measured for each of 
the 17 cuts and for several values of the coefficients across the L\- or L3-lines the amplitude 
distribution D(|^4|) and the phase gradient peak-to-peak spacing distribution D(p). We con- 
jectured that defects occur if and only if p > Psn- Indeed, we observe that in 11 out of 17 
points such conjecture is fulfilled. On the remaining 6 points the theoretical conjecture leads to 
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Figure 18. Space of (c\, C3)-coefficients for the CGLE with the Benjamin- Feir-Newell curve 
(dot-dashed), the lower bound Psn — ► °o (f Uil ) an d the stability-based estimate (dotted) for 
the transition from phase to defect chaos. The symbols refer to our numerical data for the 
appearance of defects. 



an estimation of the transition lines within a maximal error bar of 3%. The points determined 
following the conjecture are indicated as empty circles in Fig. ||. The small deviations may have 
different reasons, that we summarize below: 

(i) If fluctuations occurring during the phase chaotic dynamics are only moderate, such as 
happens near the L3 transition line or for small system sizes, more complex coherent structures 
can survive for a short time. Here we analyzed only the shortest coherent structures characterized 
by a single hump. We believe that this is sufficient to understand the main aspects of the 
dynamics of large systems. However, longer combined MAWs with more than one hump emerge 
from periodic MAWs via period doubling bifurcations. The existence of the long combined 
MAWs is limited by saddle-node bifurcations analogously to single MAWs, but these bifurcations 
occur at slightly bigger values of the parameters c\ and C3. Therefore the appearance of these 
more complicated structures can delay defect formation even if one inter-peak spacing within 
the structure is bigger than Psn of the single MAW. 

(ii) Near the L\ line the dynamical fluctuations in the phase chaotic regime are stronger than 
in the proximity of the L3 line. In this case and for sufficiently high values of the parameter 
ci we observed situations where not only the structure with the longest inter-peak spacing but 
also the neighboring structures were involved in the defect formation. 

(Hi) The assumption to consider MAWs with v = is only an approximation. If the average 
phase gradient locally (on scales P) deviates from then the saddle-node bifurcation slightly 
shifts towards smaller coefficients plf . 

As far as the (numerically) improved L3 and Li lines are concerned, we observe that both 
these lines lie to the left of the ones determined in earlier numerical studies [19]. This is due 
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to the fact that our simulations are of longer duration then those performed previously. This 
confirms the expectation that such transition lines will shift towards the Benjamin- Feir-Newell 
curve for increasing systemsize and integration times [21]. Moreover, some authors claim that 
indeed in the thermodynamic limit L\ and L3 will coincide with the Benjamin- Feir-Newell curve 
and the phase chaos regime will disappear [27]. On the basis of our simulations we cannot 
exclude such a possibility for higher space dimensions, but based on the results presented in 
this paper we conjecture that the saddle-node line for P — >oo provides a lower boundary for the 
transition from phase to defect chaos in the one-dimensional CGLE. 

Final Remarks - We have presented a systematic study of modulated amplitude waves (MAWs) 
in the complex Ginzburg-Landau equation (CGLE). These periodic coherent structures origi- 
nate from supercritical bifurcations from the homogeneous oscillation of the CGLE due to the 
Benjamin-Feir instability. The range of existence of MAWs is bounded by saddle-node bifurca- 
tions occurring for values of c\ and C3 that depend on the period P of the MAWs. Approaching 
the transition from phase to defect chaos, near-MAWs with large P occur in phase chaos, and 
defects are generated when the period of these near-MAWs becomes larger than the spatial 
period Psn of the critical nucleus. This scenario is valid for both the L\ and L3 transition. The 
divergence of Psn for coefficients in the phase-chaos regime led us to conjecture that there is a 
lower bound for the transition from phase to defect chaos. Considerations of the linear stability 
properties of MAWs in light of their tendency to increase or decrease the typical period p in 
phase chaos, has led us to a fit-free estimate of the apparent transition from phase to defect 
chaos that fits the numerical data well. 

Altogether, our study leaves little space for doubt that the transition from phase chaos to 
defect chaos in the CGLE is governed by coherent structures and their bifurcations. From 
a general viewpoint, our analysis shows that there is no collective behavior that drives the 
transition. Instead, strictly local fluctuations drive local structures beyond their saddle-node 
bifurcation and create defects. 

Outlook - We want to stress here that the extension of the analysis to MAWs with nonzero 
average phase gradients [p^ |, will be of considerable interest for experimentalists, because in 
some recent experiments concerning Rayleigh-Benard or Marangoni convection in quasi-one- 
dimensional geometries, supercritical Eckhaus instabilities of plane wave trains and the corre- 
sponding emergence of stable saturated MAWs have been observed |0|-15|. These states are 
analogous to what happens for the Id CGLE when phase chaotic solutions with v 7^ are 
considered [22,23|. 

The relevance of MAWs for two-dimensional structures is suggested by recent experimental 
evidence of MAWs observed in connection with superspiral and spiral breakup occurring in a 
Belousov-Zhabotinsky reaction [11|. Moreover, in the phase chaotic regime of the 2d CGLE the 
correspondence between long inter-peak spacings (here diameter of cells) and the strength of the 
local modulation has already been noticed numerically [27]. Additional mechanisms present in 2d 
remain to be explored. Thereby it might turn out that phase chaos exists in the thermodynamic 
limit in Id only but not in 2d as previously conjectured p7| . 
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